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, Power law is one of the the simplest forms of the relationship between different variables of a 

' system. It leads naturally to the introduction of compound parameters describing physical properties 

^SJ ' of the system. Often one of the variables of interest is the object dimension, or time. The prevalence 

of a simple power law over the entire range of dimensions or times can be helpfully interpreted as size 
or time independence of the corresponding compound physical parameter of the system. However, 
it is also often found that a simple power law only persists for some extreme values, e.g. for very 
'^s I large and/or small sizes, or very short or long times. Transitions between regimes of different 

. power law asymptotic behaviour are encountered frequently in the description of a wide variety of 

' physical systems. While asymptotic power law behaviour may often be readily predicted, e.g. on 

dimensional grounds, the evaluation of the relationship between system parameters in the transition 
range usually requires laborious detailed solution. To obviate this difficulty we introduce, on rather 
Q ' general basis, something we refer to as the merging, or 'knee' function. The function possesses 

C/3 ' sufficient flexibility to describe transitions of various sharpness. To demonstrate its usefuness, the 

merging function is applied to a variety of well-known scaling laws in the mechanics and strength 
of materials and structures. 
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I. POWER LAW MULTI-SCALING IN DEFORMATION AND FRACTURE 



^ Power laws play an important role in the study of size dependence of various physical properties. Barenblatt and 

Q Bazant '16] use functional analysis to demonstrate that, in the absence of an inherent length scale, the dependence 
^ of an arbitrary physical parameter of the system on its size must obey a power law. If two physical quantities are 
'— related by a power law, y — Cx", then the parameter C = yx~°' can be thought of as invariant in x. If x represents 
^ • system size, then C is said to be size independent, or scale invariant. This parameter is often thought of as a property 
\ of the system, and it dimension is compounded from those of the quantities x and y. 

• However, in practice it is the deviation from the simple power law descriptions that is of particular concern. The 
' terms 'size effect' or 'scale dependence' are conventionally used precisely when a particular power law no longer applies. 
That is, the parameter C, the system property, is no longer constant, and has to be thought of as a function of x. 
OO In the present context we focus the attention on the size dependence of material strength. The precise nature of 

the physical parameter y representing 'strength' depends on the context: it could be a critical value of stress, or stress 
rate, stress intensity factor, hardness, etc. We consider several instances of power multi-scaling, i.e. the transition 
between distinct power laws that persist over certain limited ranges of size. We then pose and answer the key question 
■ about efhcient description of this type of multi-scaling. 

By way of illustration, consider a plane strain crack of length 2a in an infinitely extended elastic solid subjected to 
a nominal remotely applied stress a the derivation of the Griffith P| propagation criterion begins with the following 
' expression for the total system free energy 

a : 

O. J^(a,a)=47- ^ ^ W . (1) 

• i-H It is of particular interest to note the presence of two terms in the above expression displaying the dependence 
on different powers of the crack dimension, a. The first of the two terms represents the additional surface energy 
^ , associated with the freshly created crack surfaces, proportional to the surface energy density, 7. The second term is 

■ - - ' negative and represents the reduction in the strain energy stored in the system due to the load relieving effect of the 
crack. 

The presence of the two energy terms with different size scaling behaviour ensures that the free energy of the system 
is not a monotonic function of the size. The energy rather displays an extremum, found by setting the derivative with 
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respect to size to zero, 

— = 47 ^-a^a = 0. (2) 

da ' E ^ ' 

Hence, the critical combination of crack size and applied stress for crack propagation is found to be 

2 

(T na= 2 = ^ic- (3) 



where Kjc is considered to be a material property. This compound material parameter has unusual dimension of 
MPay^, and was found to be constant for defects of sufficiently large size. 

According to the linear elastic fracture mechanics (LEFM) criterion, the scaling law for the failure stress as a 
function of crack size is given by 



'Kic 

Consequently, for specimens with cracks of length a and a\ the ratio of strengths can be written as 



(4) 



(5) 



a la 

This is the most widely known scaling law in linear elastic fracture mechanics: the failure stress decreases as the 
inverse square root of the crack dimension. This power law behaviour is represented on the bilogarithmic scale of 
failure stress vs specimen size by a straight line with gradient —1/2. 

It is immediately apparent that a paradox arises when cracks of progressively smaller size are considered. According 
to the above LEFM scaling law, the specimen strength would continue to increase without limit, as smaller and smaller 
crack sizes arc considered. This is clearly in contradiction with the observations, since even most carefully prepared 
and defect-free specimens of any material always possess finite strength. For very small defect sizes, and particularly 
when no discernible flaw can be detected, a stress-based failure criterion must be used instead: 

(7 = (To, (6) 

or, in the explicit form of a power law relationship between a and a, 

a = aoa°. (7) 

The distinction between the two extremes of behaviour within the framework of LEFM is sometimes referred to as 
the long crack and the short crack regimes, respectively. In the long crack regime the inverse square root scaling with 
crack size is expected. In the short crack regime the specimens must possess a limiting stress value associated with 
alternative failure phenomena, e.g. yield. In order to capture the strength scaling behaviour over the entire range of 
crack length, one requires a combination of the two criteria, or a multi-scaling description. 

Some simple but fundamental questions arise. 

Firstly, given that each power law applies in its respective extreme of flaw size, namely, the stress-based failure 
criterion when a << a*, and the fracture mechanics based criterion when a >> a* , how can the transition between 
the two power law regimes be described? 

Secondly, in formulating the above statement one is forced to introduce the parameter a* in order to express the 
asymptotic behaviour. What is the nature and the value of parameter o*? 

These questions arc typical of many situations in science and engineering, and arc usually given ad hoc answers 
that depend on the nature of the system being studied. In particular, the usual procedure is to develop a more refined 
and detailed model that is capable of describing the 'interaction' between the two asymptotic regimes, but invariably 
requires significantly greater effort to solve. Requirements naturally applied to such a model would certainly include 
the conditions that in the respective limiting regimes the model should converge to the simple power law descriptions. 

The approach proposed in this paper aimis to obviate as far as possible the need for developing detailed solutions 
for the transition region. Rather, we ask whether it is possible to make some simple assumptions about the nature of 
the transition and to derive a general form of the 'merging function' between the two regimes. When once such form 
is established, the merging function can be calibrated using limited amount of experimental data, and will provide a 
very simple and portable expression of the multi-scaling power law response of the system. 

In the following section we develop a form of the merging function on the basis of some fairly arbitrary assumptions. 
We do not contend at this point to claim that this is the correct, the best or the only form of such function that 
can be identified. Instead we seek to obtain a formulation that can be applied to a variety of situations in order to 
check its capability of describing the multi-scaling power law transition. We then focus our attention on illustrating 
the applicability of the proposed function in a wide variety of situations drawn from the context of deformation and 
strength of materials. 



3 



II. 'KNEE' FUNCTION DESCRIPTION OF THE TRANSITION BETWEEN TWO DIFFERENT 

POWER LAW REGIMES 

Consider a physical system described by a relationship between two parameters, x and y, that dispays the following 
asymptotic behaviour: 

x«xo, {y/yo) = [x/xoY, 
x»xq, {y/yo) ^ {x/xo f . 

The exact nature of the transition is not known. However, to make further progress we shall assume that there exists 
a non-dimensional parameter of the system, in the form of a product of powers of y/yo and x/xo, which is given 
simply by a sum of two power law functions of x/xq: 

2{y/yor = ix/xoy" + [x/x^yf". (9) 
Expressing the original parameter, y, from this form, we obtain: 



(y/yo) = {x/xoT 



1/7 



(10) 



The division by 2 within the brackets ensures that the transition between the two regimes occurs at the point x/xq — 1, 

y/yo = 1- 

We further consider a special case when one of the power law regimes corresponds to a constant, i.e. a — Q. In this 
case it is convenient to describe the scaling behaviour by the function: 

{y/y,)^[{l + {x/xoV')]"\ (11) 
One particular situation often encountered in practice is when both [3 and 7 are negative, (3 = —b, 7 = — c, leading to 



[{l + ix/xof-f" 

The above equation often describes the intraction between two different physical mechanisms that results in an increase 
of a particular property y (strength, say) over the base value, yi. It is then appropriate to write th above expression 
in the form 

{y-vi)/va = - — (13) 

[{l + {xlx^r)f'' 



leading to 



y = yi H rr- (14) 

[{i+{x/x,rf'^ 



The formulae derived in this section are rather simple, and contain a relatively small number of parameters. 
However, they are found to describe successfully a wide range of multi-scale power law relationships in the science 
of material strength. In the following sections we provide illustrations of this fact, and then discuss implications and 
draw conclusions at the end. 



III. THE TRANSITION BETWEEN STRESS AND TOUGHNESS-CONTROLLED STRENGTH 

The existence of the two regimes of strength scaling described in the introductory section, namely, the toughness and 
stress controlled strength, already became apparent in the early experiments of Irwin j^]. He conducted experiments 
on thin (0.8mm) large sheet of aluminium 7075 alloy containing central slits normal to the tensile loading direction. 
Irwin calculated the net section stress at instability, accounting for finite width of the sheet, and plotted this stress 
as a function of the crack (slit) length 2a. Irwin's results clearly demonstrated the inverse square root behaviour 
predicted by LEFM, as illustrated by the dropping portion of the curve shown in Fig. 2. However, it is also apparent 
from Irwin's results that for shorter cracks there exists a limiting value of stress that is independent of the crack 
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FIG. 1: Illustration of multi-scale power law transitions described by the 'knee' function with different sharpness parameters. 



length 2a, provided it does not exceed some threshold value, 2a < 2a\. The two failure criteria with different scaling 
behaviour need to be combined to obtain a description that is valid over the entire range of crack lengths considered. 

Fig. 3 contains the Irwin data re-plotted in the bi-logarithmic scale to highlight the power law multi-scaling be- 
haviour. The data markers appear alongside the curve that represents the merging function of the type introduced 
in the previous section. The form of the function is chosen to satisfy the asymptotic behaviour for exteme values of 
the crack length parameter a, and is given by 
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CTi (l + (2a/2ai)»)i/(2n)- 



(15) 



The best values of the other parameters appearing in this expression are given in Fig. 2 and Fig. 3. 

The two-criterion failure condition introduced here occupies an important place in fracture mechanics analysis. It 
arises naturally within the elasto-plastic fracture process zone model, that is alternatively known as cohesive zone, 
bridging zone, or strip yield model, and is common;y associated with the names of Barenblatt Q and Dugdale j^. 
This model is an important tool in the analysis of size effects in fracture. 

For our present purposes we quote the following result from the Dugdale model |^ : 



a 



exp 



%a\a 



(16) 



Here we denote by a\ the limiting (yield) stress, E* is the plane strain Young's modulus, and (5c is the critical crack 
tip opening displacement (material property). Introducing a characteristic crack length 



ai 



8(71 



(17) 



the failure criterion can be re- written in a convenient form allowing direct comparison with other expressions considered 
here, 



G 



2 

— arccos 

TT 



exp 



(18) 



We now wish to see how well the two-criterion failure behaviour can be described by the use of the merging function, 
or 'knee' function, introduced above to describe Irwin's results in Fig. 2 and Fig. 3. We use the form 
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(l + (a/a2)«)i/(2n)' 



(19) 
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FIG. 2; The Irwin et al. ^ data for the net failure stress as a function of the total slit length 2a, plotted in linear coordinates, 
together with the merging function describing the transition. 
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FIG. 3: The Irwin et al. 3 data for the net failure stress as a function of the total slit length 2a, plotted in bi-logarithmic 
coordinates, together with the merging function describing the transition. 



which is different from that of equation 115|l only in that the reference crack length ai has been replaced with 02. 
The comparison between the Dugdale solution and the 'knee' function form in shown in Fig. 4. Excellent agreement 
is found between the two descriptions, equations H18|l and (|19|l . with the following values of the parameters 

n = 1.996 ±0.052, aa/ai = 0.75 ± 0.01. (20) 

Therefore, equation (|18|) is also well approximated by the following simpler expression 

- = ^ Tu- (21) 

(1 + (24(Tia/3^i;*5c)2)^ 
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FIG. 4: The comparison between the two-criterion strength scaling according to the Dugdale model (solid line) and the 'knee' 
function approximation (dashed line). 
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FIG. 5: The experimental data of Mecholsky on the scaling of fracture zones in a ceramic described using the 'knee' 
function. 



However, equation (|19|l offers greater flexibility, e.g. in terms of describing the sharpness of the transition: in fact, we 
have already seen (Fig. 3) that the best choice of parameter n for the description of Irwin's data is n = 10, not n = 2 
that corresponds most closely to the Dugdale solution. It is also worth noting that the characteristic crack lengths in 
the two descriptions, oi for the Dugdale model and 02 for the 'knee' function, are close in magnitude and only differ 
from each other by a multiplier of the order unity. This is only to be expected from the introduction of the 'knee' 
function in the previous section, since the parameter xq in equation (|12|l has a value close to the vlaue of x at which 
the switch between two power law behaviours takes place. 

Having succeeded in matching the Dugdale curve for ductile materials using the 'knee' function, we can turn our 
attention to the case of quasi-brittle materials. Mecholsky et al. ^3 investigated the fracture surfaces of a ceramic 
(Pyroceram 9606) and identified a system of concentric annular zones emanating from the point of initiation, that 
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were identified by the terms 'mirror', 'mist', and 'hackle', to indicate the varying degree of surface roughness. The 
variation in roughness is clearly associated with the change in crack propagation speed, which in turn is likely to 
depend on the amount of excess strain energy available for crack propagation and acceleration. Since this excess 
strain energy is clearly itself a function of crack size and applied stress, it is natural to find that the radii of zones 
exhibiting certain roughness scaled with the fracture stress. Mecholsky used a simple power law relationship, 

ar" = M, (22) 

where r may refer to the radius of the 'mirror', 'mist' or 'hackle' region. For large zone sizes the above equation shows 
LEFM scaling with m = 0.5. However, plotting the failure stress a versus the radius r in bilogarithmic scale (Fig. 5) 
reveals a clear deviation from 'pure' power law scaling at small radii. The set of experimental data in question can 
once again be successfully described using the 'knee' function (Fig. 5). 



IV. FATIGUE GROWTH THRESHOLD AND THE KITAGAWA-TAKAHASHI DIAGRAM 



The transition between two power laws considered in the previous section in the context of static fracture strength 
is also observed, in virtually identical form, in the context of cyclic loading. Fatigue crack growth threshold is a term 
that refers to the stress level below which no crack growth is observed. It is found that, for nominally defect-free 
samples this threshold is expressed in terms of the stress range, Acq, while for samples containing cracks the threshold 
value of the stress intensity factor range, AKq. Kitagawa and Takahashi ^ proposed combining the two criteria on 
a single diagram, thus producing a map on which non-propagating cracks occupy the region below both curves. The 
transition crack length aq is defined by the relation 

Aao^/^Tao = AKq, i.e. ao — ( —r — -] ■ (23) 



V Ado 

El Haddad |3 noted that in the intermediate region the experimental data in fact lie some way below both curves, 
and suggested a form of the threshold stress range combining the two criteria given by 

^7r(a + ao) 

We note that this expression represents a particular case of the 'knee' function introduced here. Indeed, using the 
relationship of equation ()23f) the above formula can be rewritten as 

(25) 



(l + (a/ao))i/2- 

To illustrate the above discussion we use the fatigue threshold data compiled by Murakami and Endo IS], Fig. 6 and 
Fig. 7. Instead of the El Haddad formula H25|l . a general form (|19f) can be used with parameter n optimised to provide 
the best match to the data. It is then found that the data is better described by the formula 

^ Ao-Q ^ Actq 

" (l + ia/aoy/^) (1 + V^)' 



V. KINETIC DIAGRAMS OF FATIGUE CRACK GROWTH 



Fig. 8 shows the dependence of fatigue crack growth rate on applied stress intensity factor range, AK, for a titanium 
alloy (T^. The data has been presented on bilogarithmic scale with the crack growth rate, da/dN, chosen as the 
abscissa. This choice allows the 'knee' function to be applied for the description of the dependence, sought in the 
form 

Aif-Aif,, = -^%^^^. (27) 

where AKth denotes the fatigue threshold, AKjc corresponds to the critical stress intensity factor for fast fracture, 
a' = da/dN is the crack advance per cycle, and ag is the reference value of this parameter. Parameters b and c 
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FIG. 6: The 'knee' function description of the fatigue threshold data compilation of Murakami and Endo 
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FIG. 7: The fatigue threshold data of Fig. 6 plotted as stress intensity factor range AK versus defect size a. 



describe the slope of the diagram and the sharpness of the transition. Assuming stress intensity factor threshold of 
AKth — T.SMPay^, the following values of the parameters are found 



AKic 86 ± 7MPa^A5, ao = (6.1 ± 1.3) x 10~V, 



(28) 



corresponding to the continuous curve shown in Fig. 8. 

In the central region between the asymptotics of large or small values of the stress intensity factor range the Paris 
diagram appears to contain a region of pure power law scaling. This phenomenon is referred to as an intermediate 
asymptotic [3, and is illustrated by the straight line fit to the central region of the diagram. It is this relationship 
that became known as the Paris law. It is interesting to note that it persists over several decades of fatigue crack 
growth rates. 

Crack resistance curves, or i?-curves are often employed for the description of varying apparent strain energy 
release rate required for the onset and continuation of crack propagation. Consideration of such curves is interesting 
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FIG. 8: Paris fatigue crack growth diagram from described by the 'knee' function. The power law indicated by the straight 
line represents an intermediate asymptotic. 




FIG. 9: A typical J — R diagram (data from Miserez et al., 0) 



in the present context, since they capture the interaction between external loading and the multi-scale material 
response, beginning from microscopic tip blunting processes occurring at the immediate vicinity of the crack tip, but 
progressively involving larger regions of material in plastic deformation (and hence energy dissipation) as the load 
becomes higher, or the number of fatigue cycles is increased. The morphology and size of zones of plastic deformation 
and distributed damage play a central role in determining the scaling of observed strength. For the purposes of the 
present discussion, however, we concentrate on the discussion of the shape of crack resistance curves. 

Miserez et al. 01 presented a careful study of crack propagation i?-curves in terms of the J-integral {J — R curves) 
for several high volume fraction particulate composites (aluminium matrix containing 50% of AI2O3 or B4C particles). 
The authors identify three distinct stages of the process: 



• Initial part conventionally attributed to crack-tip blunting in unreinforced metallic alloys. Some apparent crack 
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FIG. 10: The stress-strain curve of an elastic-plastic material obeying the Ramberg-Osgood law. 

extension occurs in the composites in this study in this domain due to internal damage build-up at the crack 
tip. 

• A second stage starting close to the maximum load. This region is thought to correspond to macroscopic crack 
propagation. For finer particle sizes (10 or 5 /im) the crack often propagates in an unstable manner. 

• Fully stable crack propagation, under which a plateau value is finally attained in a region well beyond the 
domain of validity of J-controUed fracture. 

The above identification is based on mechanistic arguments. Fig. 9 shows the data of Miserez et al. re-plotted 
in bi-logarithmic scale and approximated by continuous 'knee' function curve of the type introduced here. This 
interpretation heps reveal a more detailed picture. In particular, the three distinct regions of power law scaling can be 
readily identified. Also apparent is the lack of a plateau in the extreme of large crack extensions a, and the significance 
of the range of lengths between 100 and 300 /im in the transition between different modes of material response. In the 
case considered this characteristic range of length scales corresponds to the transition from the plastic zone containing 
a small number of particles, each surrounded with locally yielded material, to the regime of extended flow with large 
numbers of particles absorbed in the plastic zone, and the material behaviour corresponding to that of structure-less 
continuum. 



In the analysis of constitutive laws of material deformation multi-scale power law behaviour plays a special role. 
In continuum deformation modelling the strain is assumed to be additive, i.e. the total strain is given by the sum 
of elastic and inelastic parts. In the case of time independent elasto-plasticity the Ramberg-Osgood law [lj| is often 
used: 



At small applied stresses the strain response is dominated by linear elasticity. However, when stress approaches a 
critical value (yield strain), the nature of the response changes, and for large stresses the strain is dominated by 
power law hardening. Ramberg-Osgood law for the description of this transition is popular both for its simplicity and 
due to the fact that it provides adequate approximation for a wide range of strains. We wish to explore the general 
background to this situation, and attempt to identify the functional form suitable for the description of the transition 
between power law asymptotics in elasto-plastic deformation. 



VI. STIFFNESS OF AN ELASTO-PLASTIC MATERIAL 




(29) 



11 



10' 



Oh 4 

o 



10^ 



10 



— Tangent modulus 
Secant modulus 



10"^ 10"^ 10"* 10"^ 10"^ 10"' 10° 



Strain, s 



FIG. 11: The tangent and secant stiffness of an elastic-plastic material. 
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FIG. 12: The reduction of the instantaneous stiffness of an elastic-plastic material with respect to its Young's modulus as a 
function of strain. 



For illustration we choose the Ramberg-Osgood description of stainless steel 316, for which the parameters were 
found to be i? = 250 ± 21GPa, A; = 436 ± 15MPa, m = 13.4 ± 1.1. The stress-strain curve for this material is 
shown in Fig. 10 in bi-logarithmic coordinates, including the markers for experimental points and the continuous 
curve representing the relationship in equation H29|l. For strains below about 0.001 linear elasticity is the dominant 
mechanism, while for strains greater than about 0.001 power law plasticity takes over as the principal mechanism 
determining total strain. 

The critical value of stress a* at which the elasto-plastic transition takes place is associated with yield within the 
framework of this model, and can be readily assessed approximately. Equating elastic part of the total strain to the 
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FIG. 13: The variation of apparent hardness of a coated system as a function of the relative indentation depth. 



plastic part the yield stress and strain values are found to be 



The tangent modulus^ or the instantaneous elasto- plastic stiffness of the material, is defined as = da/de. 
The secant modulus is defined as = a/e, and can be expressed as 

— - 1 

E ^ l + (a/cr*)™' ^ 

that is, is a particular form of the 'knee' function. It is possible to use the general form of the knee function (|12|l to 
describe the dependence of the secant modulus on the total strain. The result is shown in Fig.ll, together with the 
curve for the tangent modulus. 

In order to study the variation of tangent modulus in more detail we introduce a compound parameter describing 
the magnitude of deviation of the tangent modulus from its value in the elastic regime, E' = E-E'^. The dependence 
of modulus deviation, E' ^ on strain is shown in Fig. 12, together with the continuous curve illustrating the quality of 
the 'knee' function fit to the data. Of particular interest here is the sharpness of transition in the vicinity of yield, 
and the suitability of the 'knee' function for the description of this phenomenon. 



VII. HARDNESS OF COATED SYSTEMS 



A particular example of indentation size effect on hardness that we wish to consider here concerns the variation of 
the apparent hardness of a coated system with the indentation depth. A model for the description of this variation 
was introduced by Korsunsky et al. 4J in the form 

where H is the apparent hardness, Hf is the intrinsic hardness of the coating film, Hg is the substrate hardness, 
6 = d/t represents the relative indentation depth, where d is indentation depth and t is film thickness; and So is the 
reference value of the indentation depth. The above formula represents a particular case of the 'knee' function. 

The application of the 'knee' function description is illustrated in Fig. 13 and Fig. 14 using the vast body of inden- 
tation data collected from electroplated Ni coatings on Cu by Tuck et al Fig. 13 shows the apparent system 
hardness data in semi- logarithmic coordinates, while Fig. 14 reveals the multi-scaling power law transition by plotting 
the hardness increase over the substrate, H — Hg, against the relative indentation depth 5, in a bi-logarithmic plot. 
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FIG. 14: The 'knee' function description of the hardness increase due to the presence of the coating. 

VIII. DISCUSSION AND CONLUSIONS 

Data and analysis presented in this paper cover a broad range of experimental observations from various branches 
of the science of material strength. The coverage is not intended to be complete, as this would never be possible in 
a short paper. Phenomena not dwelt upon here, but known to obey similar laws, include time-dependent plasticity 
(creep), grain size dependence of yield stress (Hall-Petch law), etc. For a more complete analysis the reader is referred 
to the relevant sections of a fuller report . 

It is also important to note that more complex situations than those considered here may be encountered, such as 
discontinuous transitions between power law regimes, dependence on multiple parameters, etc. It is hoped, however, 
that the ideas presented here may be found helpful in developing better interpretative tools for other, more complex 
situations. 
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